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Abstract 

We present a detailed derivation of the renormalization group equations for two 
dimensional electromagnetic Coulomb gases whose charges lie on a triangular lattice 
(magnetic charges) and its dual (electric charges). The interactions between the 
charges involve both angular couplings and a new electromagnetic potential. This 
motivates the denomination of "elastic" Coulomb gas. Such elastic Coulomb gases 
arise naturally in the study of the continuous melting transition of two dimensional 
solids coupled to a substrate, either commensurate or with quenched disorder. 



1 Introduction 



The understanding of defect-mediated phase transitions in two dimensions 
relies on the renormalization group study of Coulomb gases (CG). In the 
simplest examples of the 0(2) or XY model, the criticality of the Kosterlitz- 
Thouless phase transition is described using the scalar Coulomb Gas [12]. In 
this case, the charges correspond to the integer topological charges of the XY 
vortices, which interact via the 2D Coulomb (In) potential. If we perturb the 
XY model by a p-fold symmetry breaking potential (the so called clock model), 
the previous scalar CG has to be extended : the clock potential translates 
into magnetic scalar charges [TOj. These magnetic charges mutually interact 
via the same Coulomb potential, and their coupling with electric charges is 
a Aharonov-Bohm potential ^llj. This scalar electromagnetic CG has been 
studied using the real-space renormalization techniques [TUpTj . which provide 
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the critical properties of the initial clock model. Moreover, the phase transi- 
tions of various two dimensional models, such as the Ashkin- Teller model, the 
q-state Potts model, and the 0(n) model can be studied using these scalar 
CG techniques dUni. 

An extension of the scalar (electric) Coulomb gas is required in the study of 
the continuous melting transition of a two dimensional solid|T5]. This exten- 
sion is twofold : (i) the topological charges of two dimensional dislocations 
are Burgers vectors instead of integers and (ii) the interaction between these 
vector charges consist of the usual 2D Coulomb potential [i.e In interaction), 
and an angular interaction which couples the charges to the vector join- 
ing the two defect^!]. This angular interaction spoils the conformal invariance 
of the In CG. The renormalization group study of the conformally invariant 
case was achieved in ref. [S]. Studying the melting transition in the general 
case amounts to consider the perturbation by marginal conformal (rotation) 
symmetry-breaking operators of the previous conformal fixed point. The study 
of the corresponding vector CG was performed in [T6|20] . The natural exten- 
sion of this vector CG to the electromagnetic case arises in the study of two 
dimensional melting in the presence of a translation symmetry breaking po- 
tential, e.g a. coupling to a substrate via a periodic modulation of the density, 
as in Ref. [16] . Such a general vector electromagnetic CG has never been stud- 
ied to our knowledge, and it is the purpose of the present paper to derive the 
RG equations describing its scaling behavior to lowest order. A preliminary 
study, motivated by the problem of a substrate with quenched disorder [B] was 
published some time ago, and involved a replicated VECG [7]. The present 
study provides a complete and general derivation of the RG equations valid 
for any type of substrate (periodic and/or disordered). The VECG studied 
here can be viewed as an extension to the vector/elastic case of the scalar 
electromagnetic CG [T7], and an extension to the electromagnetic case of the 
vector CG of [TOIEU] . As we will see, the elasticity manifests itself not only 
in the angular interactions of the electric/electric and magnetic/magnetic po- 
tentials, but also into the electric/magnetic interaction which is no longer a 
simpler Aharonov-Bohm potential. 

Before turning to a more precise definition of our model, let us mention the 
field theoretical approach to the CG problem. The scalar electromagnetic CG 
admits an equivalent Sine-Gordon field theoretical formulation [T9]. Its scaling 
behavior in the electric case was derived in Ref. [Ij. Extension to the electro- 
magnetic CG case were considered in p[5] (see also [5]), which included in 



This angular interaction is a manifestation of the microscopic nature of the 
dislocations, which can be viewed as additional half-hne of atoms inserted in the 
lattice[13]. A pair of dislocations of opposite Burgers vectors, which is an extra 
segment of atoms, has obviously some preferred orientation with respect to the 
initial regular lattice. 
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particular parafermionic operators[9]. This electromagnetic In — CG was ex- 
tended to consider charges in higher groupsj^, as well as relations with string 
theory models. In these generalized Toda field theories, the CG charges appear 
as root vectors of Lie algebra, and the charges of the SU{3) Toda field theory- 
can be identified with Burgers vectors of a triangular lattice. In this perspec- 
tive, our present study corresponds to an extension to the non-conformal case 
where angular interactions are included of the SU{3) study of Boyanovsky 
and Holman . 

The paper is organized as follows : in section [2], we derive the CG formulation 
of an elastic solid coupled to a substrate. We consider explicitly two impor- 
tant cases : the case of a periodic commensurate substrate, and the case of a 
random pinning substrate. This allows to define the general vector "elastic" 
CG which is the subject of this paper. In section [3l the renormalization group 
equations for this general CG are derived to order one loop, using a real space 
procedure similar in spirit to the method described in [T7j. The results are 
summarized in Section |H Finally, in section O these equations are restricted 
to the original elastic models. Due to the complexity of the present derivation 
we have deferred to a separate publication the study of these RG equations 
for the various models. 

Notations 

Throughout this paper, we use the notations Jf^ = J cPr = dxdy and 
J^= J (Pq/{2TiY. The notation r corresponds to vectors in the two dimensional 
plane, originating from either the direct or dual lattice, while boldfaces A 
denote vectors in the replica space. Vectors both in replica and two dimensional 
plane Ai^a are denoted A. The sum over repeated (real space or replica) indices 
will be assumed : 

n 

A-i^aBi^a = ^ ^ Ai aBi^a (1) 
i=l,2 a=l 

and we use the convolution notation 

[A*B]{f)= f A{f)B{f-r) (2) 

Jr' 

which for a density of charges 6(r) = J2a baS{r — ra) reduces to 

bi * Vij *bj = ba,iVij{ra - ^/3)&/3J (3) 

ij=l,2 a/3 

Unless otherwise stated, the indices i,j,k,l will correspond to real space in- 
dices i = 1,2; a, b, c, d to replica indices between 1 and n; and greek indices 
a, (3 label the charges in a collection of charges. The notation e corresponds 
to the unit vector e/\e\. 
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Fig. 1. Representation of a hexagonal lattice we will consider in this paper. The 6 
vectors ±6^=1^2,3 are the unit vectors of the original lattices (here the lattice spacing 
has been set to gq = 1), and the 6 unit vectors ±Gi=i,2,3 lies on the dual lattice. 

2 The model 

2.1 Elastic description of a pinned two dimensional crystal 
2.1.1 Two dimensional elastic energy 

In this paper, we will consider a crystal with hexagonal symetry (see Fig{T]). 
For such a lattice, the elasticity is isotropic, and the elastic energy is given by 
the harmonic hamiltonian[l3| 

Ho[u] = 2 y Uij{f)CijkiUki{r) = - J d'^f{2fiul + \ul,}j (4) 
= \j ^,n,m^AQ)uA-Q) (5) 

with Cijki = fJ'iSik^ji + 6ii6i^ + XSijSki where X, fi are Lame coefficients, and 
the tensor Uij is defined bjo Uij = ^{diUj + djUi). For later convenience, it is 
useful to define the local stress tensor cTj^ = CijkiUki = 2nuij + \5ijUkk- The 
elastic matrix ^ij{q) is given by 

^ijio) = + Cmq^Pljiq) ; cn = 2/i + A ; cgg = (6) 

^ Note that we have neglected the nonlinear component of Uij. 



where cn, cee are respectively the compression and shear moduhi, and A, n the 
Lame coefficients of the crystaL We have used the projectors Plj{q) = (jigj, 
Piji.'i) = 'it if = ^ij — Qifij- In this expression, m is a smooth displacement field, 
which corresponds to the long wavelength distortions of the original lattice. 
Within the context of elasticity, it must satisfy the condition \u{r) — u{f + 
ei)\ <^ ao, where is one of the unit vectors of the original lattice, and qq the 
lattice spacing. To go beyond this elastic description of the lattice distortions, 
one must allow for dislocations, which are the topological excitations of this 
elastic model. 



2.1.2 Two dimensional dislocations 

A two dimensional (edge) dislocation located in is characterised by its 
topological charge called the Burgers vector ba- This Burgers vector lies on the 
original lattice, and for most of our purpose, we will restrict ourselves to unit 
Burgers vectors corresponding to one of the six Cj, i = 1, ... 6. By definition, 
this Burgers vector corresponds to the increment of the displacement field 
when surrounding the dislocation : 

u{f) dl = a^b (7) 

where the contour integral circles around fa, and we choose to consider di- 
mensionless Burgers vectors b. A collection of dislocations can be described 
by the Burgers vector density 

m=T.bJ{r-fa) (8) 

a 

This density of dislocations induces a density of strain relaxed by a displace- 
ment field Ud{f), derived in appendix |Al and given by|18] 



UdA^ = ^ [Gij * bj] (^) = H E Qij{f- r„) baj (9) 
with g,,{f) = 5,,$(f) + — e,,G(r) + ^ll^^e,,H,k{f) (10) 

Cll Cll 

The potential $(r) gives the angle between the vector r and e.g the ei vector, 
G{f) corresponds to the usual {e.g lattice) Coulomb potential and Hij{f) is an 
angular potential. We regularize these potentials with a hard cut off : using 
0{\f\ — ao) = 1 if 1?^ > ao and otherwise, they are defined as 

G'(f) = (^In (^^^ ^(|r=1 -ao) ; G{f) + i^{f) = \ni^^^ e{\r\ - a^) 

(11) 

H^A^={p^-\^^^0{\f\-ao) (12) 
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where z = + iVy, c is an arbitrary constant, and we have defined for later 
convenience the logarithmic potential In the presence of dislocations, 

the displacement field splits into the above component Udir) induced by the 
dislocations themselves, and an independant smooth phonons part Uph{r) : 
u{r) = Uph{r) +Ud{r). Without any perturbation, the usual melting transition 
is studied by performing explicitly the integral over the smooth phonons field 
in the partition function. One is left with the partitition function of Coulomb 
gas with vector charges 6q,, whose scaling behavior describes the KTHNY melt- 
ing transition. However, with translation symmetry breaking perturbations, 
this usual (magnetic) Coulomb gas must be extended to a electromagnetic 
gas, as explained below. 

2.2 Breaking the translation symmetry 

In this paper, we will consider a two dimensional crystal coupled to a substrate 
modeled by a potential V{r) coupling directly to the density p(r) of the lattice. 
This coupling adds to the elastic Hamiltonian (jl]) an energy 

Hv = fp{f)V{f) (13) 

Jr 

which explicitly depends on u instead of Uij, refiecting the breaking of the 
translation symmetry. This symmetry breaking corresponds to the situation 
where the density p(r) and the potential V{f) have some harmonics in common 
corresponding to a reciprocal lattice vector G. In the following, we will consider 
either the case of a periodic potential commensurate with the lattice, or a 
random pinning potential. In both cases, we can consider that Jf^V{f) = 0. 
We decompose the lattice density as 

p(r) = Po ( 1 - dMr) + E e^'^-^^''-"^^) J + h.o.t. (14) 
where the G are reciprocal lattice vectors. Similarly, the coupling (IT^ reads : 

(15) 

where we have defined Vg = poV{r)e'^'^'^ . Upon coarse graining (or in an effec- 
tive long wavelength hamiltonian), only the reciprocal lattice vectors common 
to V{f) and p(r) will survive. In the above equation, the primed sum is on 
these common reciprocal lattice vectors corresponding to a non vanishing V^, 
which exists in the cases considered. In the following, we will restrict ourselves 
only to these vectors in common G of minimum length. They correspond to 
the most relevant perturbations near the pure melting transition. 
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2.2.1 Periodic commensurate substrate 



In the case of a periodic and commensurate substrate, we can use the symme- 
try = V*g to rewrite the second term of (fT5|) as 

^ = jX^^^^^^{G.u{r)) (16) 

G 

As previously mentionned, we will restrict ourselves to the three reciprocal 
lattice vectors Ga=i,2,3 of minimal length |Gi| arising in this sum. We will 
also use below the unit vectors G'q,=i^2,3 = Gal\Ga\- In addition, a periodic 
substrate modifies the elastic part of the energy by generating a new term 
coupling the orientations of the lattice to the substrate. Defining the local 
orientation 6{r) = ^{dxUy — dyUx), this new term can be written as 

6H, = lJ^f{r) (17) 

where 7 is a new elastic constant. This term is non zero even in the floating 
solid phase where the direct coupling f|T6l) is irrelevant, and must thus be 
included. 



Finally, we focus on the case of weak perturbations : to first order in V^, we 
can expand the cosine coupling into 



exp 



E cos(G„.^(f))] 

\ a=l,2,3 / 



1 + 



IK 



Ga 



2T 



E 



J\Gi\in{r).u{f) 



E 



rn(r}=±Gi,±G2,±G3 

m{r).m{r) 



m(r)=0,±Gi^l,2,3 



2T 



J\Gi\'rn(r).u{r) 



;i8) 



Defining a fugacity y[0,m] for the formal charges m(?^) as 



r[0,m] =y'^-'^ withy 



2T 



(19) 



we rewrite the partition function of the perturbed lattice as 



Z=l d[u,,ir-)] I d[Mrl] (n E Y[0,rnirl] 



u. 



id]. 



r fh{f) 



iph), 



(20) 



.{Ph)i 



exp (^-^ jjO^^ exp [i\Gi\ j_^rh{r). (u^'^^ + 4^^^) (f)) 
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Plugging the expression Q for u^'^\ and integrating over the gaussian dis- 
placement field u'^P'^\ we obtain three contributions to the remaining action: 



Z= n E ^[0' ^(^] exp (S[h/h] + S[h/m] + S[m/m]) (21) 

{6(r)} \ ^ / 

The dislocation interaction is given by the usual form extended to include the 
7 coupling (see appendix : 



1 



Sm = -— luf{r)a,,iu^^,^if) (22) 



2T 

2T \C66 +7 cii / 

(24) 

where the inverse Fourier transform of 

Uq) = {^Pt, + BPl^ (25) 

was determined as 

/,,(r) = ^(1 - d'^-')Uq) = 5,,^ (in + cte) + ^J^H^A^) (26) 
providing the following expressions for the coupling constants 

2 = — f ^66(cii - cge) ^ C667 \ ^ ^ / /^(/^ + A) ^ /^7 \ ^27) 
ttT \ cii C66 + 7/ ttT \ 2/i + A + 7/ 



Note that the dislocation core energy Ec in (1241) arises from the standard 
continuum approximation ( iTTl) of the lattice Coulomb interaction G(r), and 
the use of the neutrality condition Jpb{r) = 0. From now on, the core energy 
Ec will be incorporated in a fugacity for the b charges: 

Y[b,0]=y^-^ with = e""^. (28) 
The interaction between the fh charges follows from the gaussian integration 
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over u 



(ph) 



S[rn/m] = jm{q)%'m,{q) (29) 

-^^rha-fria (31) 

a 

with the couphng constants 

K ^l^iiy 1 ^ 1 \ _ T|Gi|V 1 ^ 1 \ 

The core energy Ec will incorporated from now on into the bare fugacity 
Y [0 , m] . Finally the cross coupling comes from the last term i eE] dsn]) : 

S[b/rn] = m{ro)Qtj{ra - rp)hj{r~'i3) (33) 



+ KeejkHikifa - rp) bj{fji) (34) 



with 



^5 = . ; Kg = — ■ 35 

271 cii 7 + C66 2n cii 7 + cee 

Defining the potential 

V,,{K,, K2, f) = K^6,,G{f) - K2H,,{r) (36) 

we can rewrite the above partition function as that of a Coulomb gas with 
both electric and magnetic vector charges : 

Z= Y ]]_Y[ba,ma]expS[b{ra),rn{ra)] (37) 

{fa,b{ra),rn{ra)} ° 



^ Note that, as a consequence of the hard-core regularization of the potentials (jlip . 
one can use indifferently a sum over distinct charges (X^^.^^) or not {J2a /?) ™ 
expression below. 
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with the action 



'/3,i 

(3J 



Substrate Disorder 

In the case of a substrate disorder, the potential V{r) which couples to the 
local density of atoms of the crystal is random : its distribution will be taken 
as gaussian, with variance 

V{r)V{r') = h{r-7^) (39) 

where h{r — r*) is a short range correlator and here and below — denotes 
an average over the disorder V . The two first contributions from the Fourier 
decomposition f|T5l) are 




where aij is a random stress field, arising from the long wavelength part of 
the disorder potential V{r). It induces local random compression/dilation and 
shear stress. Its correlator is parametrized as 

(yij{r)<yki{r') = 5{f-r') [(An - 2AQQ)5ij5ki + /^Q&{5ik5ji + SuSjk)] (41) 

whose bare values, derived from f|T5|) are: 

^11 = Po^i?=o ; = ; ym = plhji^cjT^ (42) 

where po is the mean density. The second part of the disorder comes from the 
first harmonic of V{r) with almost the same periodicity as the lattice, i.e. it is 
proportional to ^/y^ the amplitude of the qc^Gi component of V{r) occuring 
in (fT5|) . Since it is not invariant under a uniform shift of u it is usually called 
the pinning disorder. The random phase field in (HTj) is uniformly distributed 
over [0, 27r] and satisfies 

{eiMr)-^.'(^-^))) = 6,y6\f- r'). (43) 
The Gi, are the first reciprocal lattice vectors (of modulus Of = IGn'^/Sa^ ). 
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The average over the disorder fields and (Jij{r) is performed using the 

rephca trick introducing the rephcated field u°'{r), a = 1, ...n, and the corre- 
sponding replicated Burgers charge b°'{r). One defines: 



Z = Z» = ny4SJexp(^*i±^) (44) 

and consider the limit n = 0. We focus on the case of weak pinning disorder 
ym, and expand the exponential of the cosine coupling in (HUj) as in (ITS]) and 
perform the disorder average in (1441) : 



exp -2^^ J2 E cos (G^.u"(f) + Mr)) (45) 

\ i/=l,2,3a=l / 

n 

= l + ymT. ^ e^^^-e^-^^-^-^V^) + O(y^) (46) 

a,b=l u=l,2,3 

= ny^+ V'[0,m^e-^I'5^'S.-''-'' (47) 

m°(r) 

The replicated m charges have initially two opposite non zero components: 

= {6a,b, - Sa,b,) with bi^b2A<bi,b2<n,u = 1, 2, 3 (48) 

However, under the fusion process of the renormalization procedure, we will 
have to consider charges obtained as the sum of these initial charges. These 
general charges will be characterized by the property J2a = 0- Their bare 
fugacity, introduced in the above formula, reads 

F[0,m'^] = v/^S^™"-"^" (49) 



To introduce dislocations one can now follow the same steps as in Section 
[2XT] splitting Ua = u^^^^ + u^f^ . The average over the random stress tensor 
( 14TI) leads to the replicated elastic matrices 

0^ = 0115"" -An ; c^^ = cee^'^" - Aee ; 7"^^ = 7^'^^ - (50) 

Hence, by the same technique as in the case of the commensurate regular 
substrate, we obtain a Coulomb gas description (l38l) of the random model, 
albeit with coupling constant Ki^ ^ which are now replica matrices involving 
products and inverses of the replica elastic matrices (l50il 
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Ki/2 = ^ (c66(cii - C66)q/ ± C667(c66 + l) ^) (51a) 

^3/4 = ((C66 + 7)"' ± Cn') (51b) 

= ^4~^(^66Cn - 7(7 + cee)"^) (51c) 

ZTT 

= ^4~^((^11 ~ C66)q/ - 7(7 + Cqq)'^) (51d) 



and thus contain information both about elastic constants and longwavelength 
disorder. The only other modification is the nature of the rh charges, detailed 
above. 



2.3 Electromagnetic Coulomb gas with vector charges 



2.3.1 Definition 

To study the scaling behaviour of the two above models with and without 
disorder, it appears necessary to consider a general electromagnetic Coulomb 
gas with vector charges. In full generality, we will consider replicated charges 
6", of n components. Each component of the Burgers charges 6" lies on 
the direct lattice, while components of the charges are reciprocal lattice 
vectors. Any additional condition on the allowed charges, specific to the model 
considered, will be detailed at a later stage of the study. Our derivation of the 
renormalization equations will stick to the most general model. The partition 
function of this Coulomb gas is defined by 

Z= E n^[ba,AJexpSfc(r-;),<(r-;)] (52) 

where the sum counts each configuration of indistinguishable charges only 
once. These configurations correspond to electromagnetic charges 6^,m^, la- 
belled by the index a, both located in r^ which belongs either to a lattice 
(lattice Coulomb gas) or to the continuum plane with a hard core constraint 
(see eq. (ITT!)). These configurations satisfy a neutrality condition : 

^6;^ = ^m;^ = for each a = l,...,n (53) 
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The action of this Coulomb gas reads 

(54) 

where the interaction potentials Vij has been defined in (l36l) . and the coupling 
matrices Ki in section (12.2. ip for the commensurate potential, and in fl5T]) for 
the pinning random potential. We also define the geometrical factor 

A$ = flolCil- (55) 

where A$ = 47r/v^ for the triangular lattice, and A$ = 27r for the square 
lattice. Defining charge densities as 

6-(f) = ^6^5(f-r-;) ; m\r)=Y^mlb{r-r:) (56) 



we can express this partition function as 

(P'r 

exp ^ y 

{b''(r),m"(r)} 

witlH 



Z= exp ( / ^lnF[6"(f),m'^(r)] ) exp 5[6'^(r), m'^(f)] (57) 



This is a vector generalisation of the 2D scalar electromagnetic coulomb gas 
and of the electric vector coulomb gas which enter the standard study of 
melting. 



^ Note that the angle $ being defined up to a constant, the model is defined for 
configurations satisfying X]oSa^a-"^a ~ ^- This condition is satisfied in a bare 
model consisting of a collection of purely electric (m = 0) and purely magnetic 
(ba = 0) charges. Without this condition, a change of definition of the angle $ — > 
$ + ^0 is accompanied by a redefinition of the fugacities for composites charges : 
m] m] exp[— i0o&-''7i]- 
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2.3.2 Electromagnetic duality 



In 2D coulomb gas, the Kramers- Wannier duality corresponds to the inter- 
change of electric and magnetic charges : 6 m. In the usual scalar ECG, 
this corresponds to the interchange of strong and weak coupling regimes of 
the theory : g 1/g where g is the coupling constant of the ECG. For the 
present general VECG, this duality transformation can be inferred by by writ- 
ing explicitly the action (!58|) as 



S[b{f),m{r)] = ^ E WriK^Mrap) - ((&^r„^)(&^.r„^) - ]:{t-h% 



+ ^E[«-^^)^*(^"'"/^) 



with the convention Inspection of the above expression, and the 

relatioiT^ fjfj + f^fj- = Sy (or Hij{f-^) = —Hij(r)), shows that performing 
the simultaneous change: 



[ha, Aa) 



(59) 



and 





.K[ = 








(60) 


K2- 








= -K2 


(61) 




-^5 = 








(62) 



leaves the action unchanged. This is the duality transformation. Note that the 
symmetry by orientation change B —B (or time reversal) corresponds to 
i — >• —i. It affects only the h/fh interaction. 



3 Renormalization of the Coulomb gas 



The renormalization of this electromagnetic Coulomb gas goes along the lines 
of the Coulomb gas with scalar charges (Nienhuis) : upon increasing the real 
space cut-off ao — >■ 006*^' (corresponding to the size of the charges), we have 

^ Note also the useful relation rifj- — rjV^ = — ejj 
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to consider three different processes : (i) the simple rescahng of the parti- 
tion functions's integration measures and the Coulomb interaction, (ii) the 
screening or annihilation of charges, corresponding to the modification of the 
Coulomb interaction of distant charges by two opposite charges distant by 
less than the new cut-off aoe'^\ and (iii) the fusion of charges when two non- 
opposite charges distant by less than the new cut-off have to be considered as 
a new single charge at the new scale. We will consider successively this three 
processes. 



3. 1 Reparametrization 

Simple rescaling of the cut-off a^e'^^ into the integration measure {(Pr/a^ 
and and the Coulomb interaction (from the terms containing ln(r/ao)) results 
in the eigenvalue 

diY[h, m] = (^2 - ^ (b^.P'Kf + rh''.m!'Kf' + 2im1eijh]K^^')^ F[b, m] (63) 

3.2 Fusion of charges 

We consider the situation where two charges (bi, mi) and (b2, m2) located in 
ri and r2 are distant by less than the rescaled cutoff : < |p| < aoe*^' where 
we define p = ri — r2- The part 5*12 of the action (l58l) involving these two 
charges can be decomposed into their mutual interaction and the interaction 
with the rest of the charge configuration 5*12 = 5*1,2 + J2a^i,2 5'i,2/a • From now 
on, we will use the notation 

V^l, = V.j{Kf,K^') ; V^l^ = V,,iKf,Kf) (64) 
= + ^kjV^kiKf, O (65) 

With this notation, the mutual interaction between charges 1 and 2 reads 

+ ^ K.ag'(p16t +m^,,ag'(-p)6t.) (66) 
Similarly the interaction between this pair and another charge a is written as 

+ ^ {rnl,gt^{n - r^)hl^ + ml,gtj{r^ - nK,) + (1^2) (67) 
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The part of the partition function involving the two charges (bi,mi) and 
(b2, can be written a£3 

= E (n/ ^) 

n y[ba,Aa]l^[bi,mi]r[b2,m2]e^'''+^'^^i>2^^.2/. (gs) 

We are interested in the correction of order dl coming from this partial parti- 
tion function. To proceed, two cases must be distinguished : either the total 
charge in non zero, or bi + b2 = mi + m2 = 0. The first case corresponds to 
the fusion of charges considered below, and the second to the annihilation of 
charges (or Debye screening of the interactions), which will be considered in 
the next section. 

In the first case we have bi + b2 7^ or/and mi + m2 7^ 0. This gives af- 
ter coarse graining a non zero effective charge located in i? = {f\ + r2)/2. To 
proceed, we assume a low density for the Coulomb gas, which amounts to con- 
sider that all interdistances — rp between the remaining charges are much 
larger than Oq. This allows to perform a gradient expansion of the integrand 
exp [Si^2 + Ha^i,2 So) ■ The first non- vanishing term of this expansion is sim- 
ply the term of order for the fusion of charges. To this order, the correction 
fl68l) simply reads 

z,,2=di y: (n j^Y[h^Ma] 

(bi,mi) + (b2,ir.2)?^(0,0) 

y [bi, mi]y [b2, m2] (^J dpe^''^'^ e^-^L^ + 0{df) (69) 

where we have used the notation / dp for the integral on the unit circle dO^. 
The term fl69l) will correct the partition function over the same final config- 
uration of charges, including the new effective charge in R. To order in 
the gradient expansion, I]q^i,2 'S'i,2/q provides exactly the correct interaction 
between the new charge and the rest of the configuration. Thus the above 
partition function can be absorbed into a correction to the fugacity for non 
zero charges 

a,y[b,m]=^ J: ^ %,^,);(6.,m.)^[bi'Ai]^[b2,m2] (70) 

(bi,mi)+(b2,m2)=(b,m) 

where the numerical factor 

^(bi,mi);(b2,m2) = / ^^P Gxp (5[(bi , liii); (bs, ^2)]) (71) 
^ Note that the multiple integral should be restricted to the domain \ fa — Tp\ > oq 




16 



— * — * 

with the action S'[(bi, mi); (b2, 1112)] given by fl66|l with p = uq : 

5[(bi,mi); (b2,m2)] = - ih,)^{h,)';KfH,,{p) - {rh^)^{rn,f^KfH,,{p) 

+ ^{Jh,)t{h,)t^^p) - ^{,h,)tih,)]Kte,,H,u{p) 

(72) 

In the case K2 = K^^ = Kq = 0, the angular integration (I7T]) provides the con- 

strainst Y.a,i ((ini)i'(b2)i' + (m2)"(bi)°) = upon fusion, implying that the 

condition m^.b'^ = is preserved. Unlike the scalar case, this is not suffi- 
cient to forbid the generation of composite charges. For arbitrary K2-, K4, Kq, 
these composite charges will certainly be generated upon coarse-graining. 

3.3 Annihilation of charges : the screening 

Now we consider the situation of two opposite charges bi + b2 = mi + m2 = 0. 
The correction to the partition function coming from the configurations with 
these opposite charges still take the form of fl68|) . with the condition bi = 
— b2;mi = — m2- This condition implies that the ffist term of the gradient 
expansion, considered in fl69p . now only provides a constant term to the free 
energy, which we will neglect. To get the ffist non-trivial corrections to the 
system's thermodynamics, we have to consider this gradient expansion up to 
second order. To this purpose, we expand the action S'i_2/a in powers of p, i.e 
of aQ, with charges bi/2, mi/2 now located in R. In the present case the terms 
of order and 2 vanish as the pair 1, 2 is neutral, and we obtain 

Si,2/a = bl,p,d,V(^l^{R - f^)blj + ml^p,d,V^l^{R - r^)mlj 

+ I {ml,p,d,Gt^{R - f^)bl^ - ml,p,d,gt^{r^ - R)bl) + 0{al) (73) 

Expanding the second exponential to second order in oq, the correction fl68l) 
takes the form0 

Zi,2= E ( n /^^[b„,m„]U E F[bi,mi]r[-bi,-mi] 

{bc,ifia},Q7^1,2 \a5^1.2 / (bi,mi) 

/^/<M< .$^fi+i:«-.^/<.+^i:w..Je**' (74) 

J Of) Jao<|p|<aoe'" 0,0 ^ a a,l3 J 

Notice the ^ factor in front of the sum over (bi , mi ) , which accounts for the 
indiscernability of the charges 1 and 2. 



17 



~ ~* — # — * 

with Si^2/a given by fl73|) and Sffcij-mi] by fl66l) with b2 = — bi,m2 = — liii 



K.^?(p1&1,+K.^^'(-p)&i,) (75) 



As explained above, the first term can be neglected as it renormalizes by a 
constant the free energy. The second term vanishes by the symmetry p —p 
of the integral. Using / d^p/a^ = dl J p (where the last integral runs over the 
unit circle), the correction from (fMl) that we will focus on can be written 
exphcitly a^ 



Zi,2= fn /^^>^[ba,mj j ^^^^[(ba,!!!^); (b;3,m;3)] 

(76) 

with the (correction to the) action 



dS[(ha, ma); (hp, mp)] = 

dl^ J2 l^[bi,mi]F[-bi,-mi]y'ci2^y'ci/)e^[^i'™^]p,pi 

(bi,mi) 



(77) 



X 



red 



icd id 



■>dcjc 



^ Note that similarly to the case of the scalar Coulomb gas, the term a = /? in this 
sum generates a renormalisation of order to the fugacity y[b, lii], which will be 
neglected in the present study. 
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This correction to the action between two charges can be rewritten as 



+<Ai 



(78) 



where we define the tensors relative respectively to the integration over p and 
R : 



dl 



(bi,mi) 



dl -* r - ^ 

[M2]:i, >^[bi, ifii]i^[-bi, -ifii] j dp e'^'^'^^^psp.bl.ml, (79b) 



(bi.mi) 



dl -* f - ^ 

[M3]Tt,k-j i: Y[h,,A^]Y[-h,,-l^,]J dpe'^'^^^^^p.ptml,^^^^^ (79c) 



(bi,mi) 



[II 



(l,3)-\ah;cd 



J s,ij;t,kl \ 



(l)]ab;cd 
2 Js,ij;t,Wl' a 



- r-» = j d'R dsV^l^iR - fo)dtgi^{R - fp) 



(80a) 
(80b) 
(80c) 



where all the above expressions are symetric in a, 13. We have used that all 
derivatives of the potentials V and Q are odd. 



19 



To proceed, we thus have to (i) perform the integral over i?,, i.e calculate 
explicitly the tensors /i,2,3 (h) perform the integral over p, i.e calculate ex- 
plicitly the tensors Mi^2,3, and finally (iii) contract all the tensors in flTHj) . 
If this final contraction can be cast into contributions to the initial poten- 
tial V^^^f^yQIji^f), this will prove the renormalizability of the present vector 
Coulomb gaz to one loop. 



3. 3. 1 Integration over R 

We focus on the tensors /i,2,3, which are all integral of double products of 
gradients of V,Q. These integrations are conveniently done in Fourier space, 
and we start by obtaining Fourier representation of these potential's gradients 
: with the definition of the projectors -Pjj(g) = qiqj and P^j{q) = Sij — qiqj = 
^ik^jiPkiiq) ^ obtain, from the definition 



the expression or its gradient 



27T 



(81) 



d^q Aq.rq 



(27r)2 g2 



d^q Aq.rqsqtqi 



{27ry 



q- 



f82) 



Similarly, using the equality fl65l) the second gradient reads 



2tx 
-2i7r 



2-K 



-2iDf 



ab 



5"" estCijkiil, 0) + emjSstQ^kiiK^, Kq 



d^q Ag.fqtqkqi 



ijkl,st I (^27r)2 " g4 



(83) 

d'^q Jq.rqt pL 
n2 kl 



{2n) 



^5) 



where we have defined 



qUK^, K2) = (i^i - K^rd.kSji + (i^i + K2ri6,,6M - 6uSjk) 



Dtjki,st = -^S''''est[5ik5ji + 5ij5ki - SuSjk] (87) 
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With this representation, the integrals Ii,2,3 are expressed as 

(88) 

{I^^rSikii^- - ^P) = -^^'CtlniK^, K^)Dl%,,^Q^n,,su{ro. - rp) (89) 

[-^3]s,jj;t,fc/(^a ~ ^p) = ~^'^'^Dijmn,su^'klpq,tvQmnpq,uv{'i^a — ^p) (90) 



where we have defined the integral 

^klmnsATa Tfs)- J a KJ J 6 



{2nyj (27r)2 q^{q'y 



dsdtdkdidmdn J {2^^ q 
Using the Schwinger representation, the last integral yields 



^2 



(27r)2g8 Qj {2tiY Jo 12 Jo 2-ku 

The differenciation of the gaussian up to order 6 is now straigthforward : 



dsdtdkdidmd, 



-A^6stSkiSrnn + c.p. (ISterms) 

+ A'^ r^rtSkiSmn + c.p. (45terms) 
- rsTtTkridmn + c.p.(15terms) 



+ A rsnrkrirmTn 



e-^'' (93) 



where c.p. means circular permutation of the indices (the number of corre- 
sponding permutated terms is indicated). Finally, using 
1), we find : 

1 / \ 1 

12 X Qstklmn{r) =——El --Ji: ^sthl^mn + 7^ rsTtdkl^mn- 

loTT Y 4L"'y 87r 

- "T- fsrtfkriSmn + - TsTtrkrirmrn + (c.p.) (94) 



In this expression L stands for an IR cut-off. We will use the following asymp- 
totic limit for the exponential integral[6] : Ei{—x) ~ 7 -|- ln(x) in the limit 
X — 0. The expressions fl86ll87f94l) . together with the contractions formula 
fl88f89f90p constitute our final explicit expressions for the integrals /i,2,3- 
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3.3.2 Integration over p 



The invariance under 27r/3 rotations of the integrals in Mi 2,3, defined in eq. 
( 17^ . ensures that these tensors are isotropic, provided that the fugacity of 
a vector charge y[b,m] is constant under any rotation of the charge b, m 



(in particular, this implies y[b,m] = Y 
decompose the tensors Mi^2,3 according tc3 



-b, — m]). Using this isotropy, we 



[MJ,:hk = dl[{TZ-Tl^)T,,,, + VZTs,,^) ■ w = l,?> (95) 

AM2]Tt,k = dl ((r- - f-) Ust,k + n^Ust,k) (96) 

and where we used the definitions of the symetric and antisymetric tensors 

Tst,ik = ^st^ik ', Tst,ik = SsAk + Sti^sk (97) 

Ust,ik = ^st^ik Ust,ik = ^sk^it + ^tk^is (98) 

By using (note the unusual definition of the trace) : 

Ti{AB) = Ast,ikBst,ik, (99) 

st.ik 

Tt{T^) = 4, Tr(Tf ) = 4, Tr(f f ) = 12, (100) 

Tr(f/2) = 4, Tr(f/2) = 12, Tr(f/f7) = 4, (101) 



we obtain the formal expression for the coefficients V^, T 



ac 



c?/r^^ = ^Tr(TMr) ; w = l,3 (102a) 

c?/f- = -lTr(TM-) + ^Tr(TMr) ; w = 1,3 (102b) 

dlT^" = ^Tt{UM^^) (102c) 

rf/f^^ = --Tt{UM^') + -Ti{UM^^) (102d) 



^ The tensor U and U arises as can be seen e.g. by expanding the definition of 
M^^) to first order in Kq which yields tensors of the form (in the case m.b = 0) 

^ j dppsPthimk (^{p.m){p.b^) - ^^-b^^ 
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with 



Tr(TMf'=) = - Y Y^[h,in\{b\b') / dp e^^^'™] 

(b,m) 



Tr(fMf' 



Tr([/M 



rf/ ^ Y^[h,m] J dp e^^^^'^\p.b''){p.b') 

(b,m) 
(b,m) 



(b,m) 

r// ^ r - ^ 

TiiTM^') = - Y^[h,m\{rff.m') / dp e^[^'"l 

(b,m) 

Tr (TMg"^) = rf/ 5] r'[b,m] j dp e^[^''^](p.m")(p.m") 

(b,m) 



(103a) 
(103b) 

(103c) 
(103d) 
(103e) 
(103f) 



Note the following useful relations : 



Tr (TMf^) - Tr (TMf' 



Tr (TMg"^) - Tr(TM 



Tr(f/Mf ) - Tr(f/M, 



r) 



9 



d 



ac 



( \ 

dl Y y%rn] J dpe^I^'™] 

V (b,m) J 



dl Y y^ihm] dpe 



(b,m) 



d 



dl 



J[b,m] 



Y Y%rh] dpe 



V (b,m) 



oS[b,m] 



(104) 
(105) 
(106) 



3.3.3 Final contraction of tensors 



With the above expressions for the M and I tensors, we can now explicitly per- 
form the contractions of eq. fl75]) . This tedious task is performed using math- 
ematica. We find that the result can be cast in the same form as the original 
interaction with changes dKi in the couplings : this proves the renormahzabil- 
ity of the model to order Y"^. Additional constants are produced which correct 
fugacities to cubic order in Y. The result of these contractions is presented in 
the following section. 
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4 Resulting RG equations for the general model 



In this Section we collect and analyze the RG equations for the fugacity vari- 
ables y [b, m] and the matrices i^j, i = 1, ..6, which parameterize the general 
VECG model defined by the action fl55]) . 

4-1 Scaling equations for the fugacities 

The equations fl63|70l) provide the full equations for the fugacities : 

diY[h, A] = (2-^ [W.l^Kf + m\m^Kf + 2im1tijh]Kf)^ F[b, m] 

+ E ^ %.mO;(6..m.)^[bi,mi]r[b,,m2] (107) 

(bi,mi)+(b2,m2)={b,m) 

where the numerical factor mi) (b2 m2) defined in eqs. ( I7T]) and (!72|) 

^(bi,mi);(b2,m2) = / ^^P exp (5[(bi , liii); (bs, ms)]) (108) 
with the action S'[(bi, mi); (b2, m2)] given by (!66|) with p = Oq : 

5[(bi, mi); (b2, m2)] = - Kthl.H{p).h\ - KtAl.H{p).va\ 

+ z(m^b^ + m^.b?)^$(p) 

- iKt (m?.i7(p).bf + m^.iJ(p).bf ) (109) 

The evaluation of the coefficients A^g^ mi) (b2 m2) model dependent. For the 
models considered here, it will be performed in subsequent publication. 

4-2 Scaling equations for the couplings matrices 

The RG equations for the coupling constants Ki are obtained by performing 
the tensors contractions of eq. fl78|) . The resulting expression is displayed in 
the appendix [Bl Here we show that their structure can be further simplified 
by introducing the new couplings pi defined by: 





= 271 {Ki + K2) 


; V2 


= 271 {Ki 


-K2) 


(110) 


P3 


= 27r {K, + K4) 


; Pi 


= 271 {Ks 




(111) 


P5 


= 27r (7^5 + K^) 


; p& 


= 27T {K, 




(112) 
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In the general case the pi (and the Fj and Fj below) are commuting replica 
matrices. Quite remarkably, the 6 scaling equations (IB.ip decouple into two 
independent set of 3 equations for the groups Pi,P4,P6, and P2,P3,P5- 

dipi = - Tipj + TipI + 2T2P1P6 + 2f 2Pi(A0 + pe) 

+ Ts{\l+pl) + Ts{\^+Pef (113a) 
dip4 = + Ti{XI+pI) - fi(A<^ -pe)^ + 2F2P4P6 - 2f2(A<^ -p6)P4 

-TspI-TspI (113b) 
dipe = - TiPiPq + f ipi(p6 - A,^) + (F2 + f2)(-Aj +^5 - P1P4) 

- ^sP^Pg - f 3P4(A,/, + Pg) (113c) 

and 

diP2 = - F1P2 - ^ipI + 2F2P5P2 - 2f 2(A<^ + P5)p2 

+ T,{Xl + pl)-ts{\^ + p,f (114a) 
dip-, = + Fi(Aj + pI) + f i(Aj - p,f + 2F2P3P5 + 2f 2(A<^ - p,)ps 

-TspI + TspI (114b) 

dlP5 = - ^lP2P5 + riP2(A</, -ps) + (r2 - f2)(-Aj +P5 -P2P3) 

- TsPsPs + f 3(A0 + P5)p3. (114c) 

where the Fj and Fj were defined in fll02[ 11031) . Their flow equation can be 
deduced from the fugacity RG equation given in the previous section. 

In addition these equations possess remarkable symmetries. The following 
transformation: 

pi ^ P2,P3 ^ P4,P5 ^ P6 (115a) 

F, ->F,;fi^-fi,^ = l,...3. (115b) 

exchanges these two groups. In terms of the Coulomb gas couplings, it corre- 
sponds to K2 —>■ —K2] —K^; Kq — > —Kq. It can be viewed formally as a 
71 / 2 charge rotation (6, m) (6-*-, m^) in the original action. This means that 
a model where the signs of K2, K^^, Kq are simultaneously changed is the same 
(up to an immaterial global rotation) with the same fugacities. 

The second symmetry is the previously discussed electromagnetic duality. It 
operates inside each of these groups, i.e the RG equations are invariant under: 



p'i=Pi ; p'i=Pi ; p'& = -p& (116a) 
p'2=p-i ; P3 = P2 ; p^ = -P5 (116b) 

F; = F3; f ; = -f 3; F^ = -F2; f '2 = -f 2; = Fi; f', = -fi (116c) 



25 



5 Resulting RG equations for the Elastic Models 



We now focus on the models defined at the beginning of the paper, i.e. an 
elastic lattice with dislocations in presence of a substrate, which can include a 
periodic modulation and/or a substrate with quenched disorder. At the bare 
level these models do not span the whole space of the six Ki (considered in 
the previous Section) but only a "3 dimensional" subspace of Coulomb gases 
(called below the "elastic sub-manifold"). Indeed these models correspond to 
the same definitions flSTj) of the couplings constants Ki (resp. replica matrices) 
in terms of the elastic constants (resp. matrices) cn, cee, 7- We find, and this is 
one of the main results of the paper, that this sub-manifold, i.e. the definitions 
fl3T]) . is preserved by the RG flow. We emphasize that this property is far from 
obvious, and cannot be easily inferred from the structure of the RG equations 
(IB.ip of the full Coulomb gas, without any knowledge of the definitions (15T|) . 



5.1 Stable elastic sub-manifold 



Let us start by expressing the coupling constants/matrices Pi in terms of the 
elastic constants/matrices. The constants from the first group read 

Pi = 2n{K^ + K2) = ^C66(cn - C66)cn\ (117) 

P4 = 27c{K, - K,) = r|Gi| ViS (118) 

P6 = 27r(/^5 - K,) = aolGil (2^ - l) . (119) 

\ Cll / 

Note that these 3 constants depend only on cii,cqg, and not on 7. These 
equations can be inverted into 

cn = T\G,\'p^' ; cge = ■ (120) 



2ao - pe 



We recall that = aolGil 



The scaling of 7 (together with cqq) is described by the second group of cou- 
plings : 

4^,2 

P2 = 2n{Ki - K2) = ^C667(C66 + 7)"', (121) 
P3 = 27r{K3 + K,) = T|Gi|2(c66 + 7)"', (122) 
P5 = MK5 + Ke) = ao|Gi|(c66 - 7)(c66 + 7)"', (123) 
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which are inverted into 



C66 = -7^ , ; 7 = -7^ • (124 

2ao P3 2ao Ps 



From these considerations we find the equations defining the "elastic sub- 
manifold" . 

- Pe = Pm- (125a) 

XI-pI=P2P3. (125b) 

(A</, + Ps) (A</, - pe) = P1P3, (125c) 

The last relation is obtained by equating the relation (11241) with (I120p . The 
second is nothing but the first, after the it/2 rotation symmetry {b,m) 
{b^,rh-^). These equation also imply: 

(A</, + P6)(A</, -ps) = ; {P5-Pfif = {Pl-P2){P3-P'i) (126) 



It is now simple to check that the "elastic manifold" (11251) is preserved by the 
RG. For the two first conditions it is straightforward, and for the third one 
can show and use that: 

PidiPs + Psdipi + (A^ + P5)diPG - {\^ - P6)diP5 = (127) 



We can now write the RG equations restricted to this subspace. Using the 
above expressions of the pi in terms of the elastic matrices, we obtain the 
main result of the paper: 

di{cn - C66) = - (Fi - 2fi) ^ (cn - c^ef (128a) 



+ (F2 - 2f 2) 2ao|Gi| (cn - c,,) + (F3 + 2f 3) ^T\Gi 



2 



dicee = - Ti^cl, - F22ao|Gi|c66 + Ts^— (128b) 

da = - (ri + 2f,) |I7^ 

+ (F2 + 2f2) 2ao|Gi|7 + (F3 - 2f3) ^T\G,\^. (128c) 

where the Fj and Fj were defined in (I102[ 11031) . Their explicit calculation and 
analysis of the resulting equations, in the specific models, go well beyond this 
paper. 
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5.2 symmetries 



Let us comment the symmetries of these equations. They are invariant under 
the transformation: 

c'n - ^6 = 7 (129) 
Cqq = Cm (130) 

i = Cii - C66. (131) 

which, as noted above, results from invariance under a 7r/2 charge rotation 
(11151) . It means that if (cii(/), C66(0) 7(0? m]) is a solution of the RG flow, 
then (c']^^(/), Cgg(/), 7'(/), y/[b, m]) is also a solution. The self-adjoint manifold 
corresponds to K2 = = = 0, i.e. 7 = cn — cee which corresponds 
to an isotropic elastic energy and interaction between charges. It is a family 
of conformally invariant VECG. Examples have been studied in [13] (electric 
case) and in [3]. 

Similarly, the electromagnetic duality (11161) is written as 

T2|Gi|2 

c'n - = l/(cii - cee) (132) 

^6= .2 1/^66 (133) 



4a, 



y ^ ^ i/r (134, 

It means that if (cii(/), C66(/), 7(0) ^«[b, m]) is a solution of the RG flow, then 
(c']^;^(/), Cgg(/), 7'(/), F;[m-'-, b"*-]) is also a solution. Hence there is a self-dual 
submanifold invariant by the flow, defined by: 

cu ~ Cm = Cm = 1 = T\Gi\/{2ao) (135) 
ll[b,m] =y5[m^,b^] (136) 

In the space of elastic constants this self-dual point forms a " line" as T varies. 
This manifold is clearly included in the "conformal submanifold" cn — cee = 7 
defined above (it obeys Ki = K3, K2 = —K4 = 0, K5 = Kq = ). 



6 Conclusion 



To conclude, we have shown how to derive the RG equations of pinned two di- 
mensional defective solids from generalized "elastic" electromagnetic Coulomb 
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gases with vector charges, defined in fl52|54l) . These RG equations were ob- 
tained to lowest order in the charge fugacity, and displayed in full generality 
in appendix [HI They involve, in addition to charge fugacities, six elastic coeffi- 
cients (or replica matrices in the disordered case). We found that they decou- 
ple in two independent sets of simpler equations (11131) and (11141) which obey 
two additional symmetry relations. We found that these general equations ex- 
hibit a restriction to only three scaling elastic coefficients, corresponding to 
the initial pinned elastic models, which we showed to be preserved under the 
RG fiow. This provides our final result : eqs. (11281) which is still sufficiently 
general to include all known cases, e.g. the scalar electromagnetic Coulomb 
gas[T7], the scalar vector Coulomb gas describing the melting transition of 2D 
elastic solids[T5], together with various extensions, e.g. the melting transition 
of pinned 2D solids [7J. Their detailed analysis is the subject of a separate 
publication. 
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29 



A Two dimensional dislocations 



In this appendix, we derive the displacement field corresponding to a finite 
density of 2D edge dislocations (and of disclinations) in the presence of a 
coupling to a substrate. We present it here for sake of completeness, and to 
clarify the notations used in this paper. 

We consider a 2D isotropic elastic lattice coupled to a periodic substrate ac- 
cording to if = I / (Pr 2fiujj + \u\j^ + '-^{eijdiUjY . Without dislocations, the 
phonon displacement field u is single valued and satisfies eijdidjU = 0. Using 
this property, we can show that (apart from boundary terms) 




which implies that the coupling constant to the subtrate 7 can be incorporated 
in new Lame coefficients A = A — 27 and /i = /i + 7 and thus is not a new 
independant elastic constant of the lattice : 

£f 2ixul + Xul, + lO^ = \j d'r 2(/x + ^)ul + (A - 27)^^ (A.l) 

This transformation can also be written as Cn — > Cn = Cn, Cqq Cqq = C66+7. 
As we will see, the appearence of dislocations breaks this symmetry. 

The local equilibrium condition for the hamiltonian flA.ip H = ^Ui * * uj 

reads 

dH 

— = ^ Mij * Uj = 2fi djUij + A diUkk + jejiemndjdmU„ = (A. 2) 

Only for non singular fields does the matrix Mij reduce to : Mjj(q) = g^[(2/i + 
A)Pj^ + fiP^] where we have use the modified Lame coefficients introduced 
above. 

Since dislocations correspond to topological singularities of the lattice, they 
induce mult i- valued displacement fields Hence if we want to formulate the 
problem of the determination of their displacement field as a classical elasticity 
problem, we need to split the displacement field Ui into a multi-valued part 
and a smooth component Ui : Ui = +Ui. The field m| provides the necessary 
multi-valueness : = bi* where b(?^ = J2a ^(r ~ r„)bQ, is the dislocations 
density. Using 9j$ = —eikdkG, we find the contribution of uf to (lA.2p : 
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1 



^0 



— {bjeik + biejk)dkG (A. 3) 



// = ~Mij * uj = -2fi dju'j - A diulk - 'yejie„,ndjdmK (A.4) 
= ^b, [{jl - 2^)e,kd,dk + (A + 2^)e,kdA) G (A.5) 



We are now facing a classical elasticity problem consisting of finding the re- 
sponse of an isotropic 2D lattice under a local force // : Mj^ *Uk = ff, which 
can be inverted in Fourier transform as 

^,(q) = M^^(q)/;(q) 

with 

^ g2 2/2 + A ^ ^ ^ 

and 

// = - 2j)e,,P^, + (A + 27)e,,/^^,). 

We end up with a displacement field u given by 

Using the approximate Fourier transform 

and the relation ejkTiTk = ^ikrjfk + Uj, we obtain the result : 

27r V /i(2/i + A) ' ^ ^ /i(2/i + A) ' / ^ ^ 

Thus the total displacement field due to a density of dislocations b is = 
^Gij * bj where 

C11C66 CiiCge 

(A.8) 

This expression reduces to the known formula |8j without any coupling to the 
substrate 7 = 0. We also realize that in the presence of dislocations, this 
coupling 7 can no longer be incorporated into renormalized elastic coupling : 
its corresponds to a third independent constant. 

To obtain the effective interaction between the dislocations, we first express the 
strain tensor corresponding to a collection of dislocations : Mij(q) = u^j{q) + 



31 



u 



'q) where, using (1A.3P 



(A.9) 



and from (IA.6I1 



. 6, /2A + 47 pL^ ^-27 



(A.IO) 



Now plugging this strain tensor into the elastic energy and using 



Mii(q)Mij(-q) = —^bi{q)bj(-q) 



4g2 



«*iiqj«iii-qj 



2g2 
1 



2A + 47\ 

Ukk{q}Ukk{-q} = ^\^- 2il+~x ) ^'^^^^^^'^^^^^i 



PL 



we obtain the desired result : 



,2q 
1 
2^ 



:&i(q)&j(-q) 



C66 2/i+A 

„L , 4c66(cii - Cm) + ^' 



— P-^ + 



Cll 



(A.ll) 
(A.12) 



Note that another method to obtain this interaction, incorporating in partic- 
ular the contribution of disclinations, is to use the so called Airy functions [8|. 
However it does not provide the displacement field, necessary in the present 



case. 



B Renormalization Group Equations for the Full Model 



In this appendix, we present the RG equations for the full VECG. In these 
expression, the couplings Ki, Fj and Fj are commuting replica matrices. 
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diKf = --2nTi [Kl + Kl) + 47rf li^ifsTs 

+ 47rr2 {K^K^ - K^K^) - At^T^ (^K^Ke - K2K5 - j 

+ 27rr3 (kI + Kl + (^^) ' j - 47rf 3 {k,^^K,K^ (B.la) 
diKf = -47rriXiX2 + 27rf 1 (xf - ^2^) 

+ 47rr2 (/Tsi^s - i^ii^e) + 47rf2 ( ^^1/^5 - i^2i^6 + ^ 



{2n)J 

AttTsK.K, + 27rr3 I + + 2K,^^ + (^^j j (B.lb) 



+ 47rr2 (iTgi^s + K^K^) - 47rf2 ( 7^47^5 + K^K^ - 



{2n)J 

- 27rr3 (X| + X|) + 47rf3X3i^4 (B.lc) 

+ 47rr2 {K^K, + K.Kq) - 47rf 2 ^K^K, - K^-^ + K^Kq^ 

- AnVsKsK^ + 2nrs (x| - X|) (B.ld) 



diKf = 2ttT, {K^Kq - K,K,) + 27rf 1 (-^K^ - K,Ke + K2K, 



{2ny 



\ \2 



+ 27rr2 I - 1 77T^ 1 +Ki + K^- K^Ka + i^2i^4 



- 27rr3 (i^3i^5 + K^Ke) + 27rr3 ( + K^^K^ + K^K^ ) (B.le) 



(27r) 



+ 27rr2 {K2K^ - KiK^ + 2K5KQ) 
A^ 
,(2^ 



27rr2 I - I TTTT 1 + + - ^1^3 + i^2i^4 



27rr3 (K^K^ + KsKq) + 27rt, [k,-^ + KsK, + K^K^^ (B.lf) 
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